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Abstract 

We study Pontecorvo neutrino-antineutrino oscillations both in vacuum and in 
matter within a field theoretic approach, showing that this phenomenon can occur 
only if neutrinos have a Dirac-Majorana mass term. We find that matter effects 
suppress these oscillations and cannot explain the solar neutrino problem. On the 
contrary, a vacuum neutrino-antineutrino oscillations solution to this problem exists. 
We analyze this solution and available data from laboratory experiments giving 
stringent limits on and Majorana masses. 
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1 Introduction 



One of the most striking problems in Elementary Particle Physics is that of neutrino 
mass. In the Standard Model by Glashow, Weinberg and Salam neutrino are massless, 
but there is no physical explanation for this fact. On the contrary, in the framework of 
Grand Unified Theories massive neutrinos appear to be more natural and, moreover, 
when quantum gravity effects are taken into account, a non vanishing neutrino mass term 
naturally arises [§. 

Furthermore, massive neutrinos are also required in Cosmology and Astrophisics to 
explain the large scale structure and hydrogen ionization in the universe (see Q and 
ref. therein) and peculiar velocities of pulsar |^. From the experimental point of view, 
indications of a non-vanishing neutrino mass came from the so called "solar neutrino 
problem" [^] and atmospheric neutrino flux anomaly 0]. 

Actually, for example, the most simple and natural explanation for the solar z/g flux 
deficit is given in term of neutrino flavour oscillations in the solar medium f^, according 
to the MSW mechanism 0, occurring only if neutrinos are massive. 

From a phenomenological point of view, the problem of neutrino mass is strictly re- 
lated to lepton number conservation (both individual ( Lg, L^, ) and total (L) ), and 
the Dirac or Majorana nature of neutrinos themselves. In fact, if neutrinos are mass- 
less, the conservation of all lepton numbers (Le, L^, L^, L) is allowed, according to the 
Glashow- Weinberg-Salam theory, due to the invariance of the electroweak lagrangian for 
un arbitrary (global) phase transformation of the matter field. Instead, if neutrino have a 
non-vanishing mass and are of the Dirac type (particle states different from the antipar- 
ticle ones), so that to the e.w. lagrangian the following mass term is added : 

= -J2^i'rMI:',uil + h.c. (1) 

1,1' 

with /, /' = e,^,r and , in general, a complex non diagonal 3x3 matrix, individual 
lepton number is no longer conserved (L^ is not invariant for i>i e^°''h'i), while the total 
lepton number is still conserved (L^ is invariant for Ui e^^vi). However, if neutrinos 
are massive Majorana particles (particle states coincide whith antiparticles ones), so that 
the mass term is 

C = -IJ:^'rM!^uil + h.c. (2) 
^ 1,1' 

with M^, in general, a non symmetric 3x3 matrix, no lepton number is conserved ( 
is not invariant for every (global) phase transformation). Note that, while the term (H) 
does not require additional neutrino states besides those present in the Standard Model, 
the term (|l|) involves right-handed neutrino states (and related antiparticles) not present 
in tha S.M. but predicted in many GUTs (see for example 0). 

The phenomenological consequences of eq. (|l|) or (0) are very intriguing. 

For istance, if and are non diagonal (in analogy to what happens with the 
quark mass matrix), neutrino flavour oscillations are predicted fl^. But, furthermore. 



many other processes involving charged lepton, which violate lepton number conservation. 
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are allowed. For example, with Dirac neutrinos, the decays /i — 67, ;U — 3e, r — ctt" 
and the conversion /z-e in presence of nuclei, such as fi~ +Ti e" + Ti, can be realized, 
while with Majorana neutrinos neutrinoless double beta decay (Z,A) — > (Z+2,A) + 2e~ 
and reaction /i~ + Ti ^ e"*" + Ca can also occur. 

More in general, one can consider a lagrangian mass term whith both Dirac and 
Majorana terms 

- C^'' = J2v,n Ml?, + \ Y.V\,^ Ml, + \ Ml + h.c. (3) 

1,1' ^ 1,1' ^ 1,1' 

involving ul and i/'^, as well as z/^j and z/£ predicted in many GUTs. This scenario, from a 
theoretical point of view, allows to give a very small mass to neutrinos in a very natural 
way through the so-called "see-saw" mechanism |TT|. The mass eigenstates coming from 
(0) are in general Majorana states, and the phenomenology previous described in the 
discussion of eq. (H) applies in this case as well. However, as we will show in this paper, 
another interesting phenomenon (also violating the total lepton number) is now possible: 
neutrino-antineutrino oscillations. Using field-theoretic methods, in section II we study 
neutrino-antineutrino oscillations in vacuum, rederiving the formula for the survival prob- 



ability first introduced by Pontecorvo [|12|. In section III we shall extend the results found 
to keep into account matter effects for neutrinos propagating in a medium. In section IV 
we analyze the available experimental data and give limits on the u^, Majorana mass. 
Finally the conclusions. 



2 Neutrino-antineutrino oscillations in vacuum 

We add the Dirac-Majorana mass term (0) for neutrinos to the standard electroweak 
lagrangian. Since here we are interested only to neutrino-antineutrino oscillations, let us 
assume that the mass matrices appearing in (^) are diagonal and concentrate our attention 
on only one flavour at a time. So the propagation of a neutrino in vacuum is described 
by the lagrangian 

C = V /cu + v'-^ }iv'^ — -^"^D {vv ^ z/'^ I/*" ^ — ^niM {vi'^ + v*-^ (4) 

where /c^ = (ct;,k) is the neutrino 4-momentum and m£,,mM are respectively the Dirac 
and Majorana mass for the given flavour. The neutrino equations of motion are then: 

]kv — -rriDV — -uimt^ = (5) 



2 2 

/cu^ - -niDi'^ - -muv = (6) 



equivalent to 

{/: - m+) (u + u^) =0 (7) 



- m_) [-U + u^) = 
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where ^ 

m± = - {mo ± niM) (9) 

From (j^) it is easy to recognize that the mass eigenstates are Majorana states, {ly + u'") 
and {—u + z/'^), with C = +1 and C = — 1 respectively. 

Let us now adopt the chiral Weyl base for the Dirac gamma matrices and define: 

We introduce the hehcity eigenstates 

a-k 



k 



0A = A(/>A (11) 



with A = ±1, and write 



l^R = X(l)x = Z(j)x (12) 

Ul = ycpx = w(j)x (13) 



the equations of motion (H), (S) imply: 



where 



UJ + Xk 


-m+ ^ 




— m+ 


UJ — Xk J 


\ n2 


u + Xk 


— m„ ' 
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— m„ 


UJ — Xk I 




y + w 




Hi 


X + z 




n2 



(14) 
(15) 



Til = y + w Hi = —y + w (16) 

n2 = X + z n2 = — X + z (17) 



From (|Tj) (for ([T5|) it suffices the replacements m+ m_ and ni, n2 — > ni, 71,2) one gets: 

(c<;2 - P - ml) rii = (18) 
w + AA; 

?T.2 = rii (19) 

so that the dispersion relations (in vacuum) for the C = ±1 Majorana states are 

uj^ - = ml (20) 

The time evolution is described by: 

|ni.2(t)> = e-^^* |ni,2(0) > (21) 

\rH:2{t)> = e"'^* 1^2(0) > (22) 
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In the ultrarelativistic limit 



2 

mi 



a;i ^ + ^ (23) 



2 

m 



Ur, k + ^ (24) 



From (|2T|) . ( P2|) and ([T6|) we may obtain, the time evolution of the left-handed weak- 



interacting neutrino: 

\y{t) > = i (e-'"^* + e"''^^*) | ?y(0) > + \ (e"'"'* - e"''^^*) | ti;(0) > (25) 



2 

If at t = we produce (through a weak-interaction process) a left-handed neutrino, at later 
times there is a non-vanishing probability of detecting (again through a weak-interaction 
process) the corresponding antiparticle, i.e. a right-handed antineutrino. The survival 
probability of the initial ul, P{^l ^l) = \ < ?/(0) | y(t) > p, is given by 

P{ul i^l) = ^ (1 + cos(c^i - uri)t) = 

^ 1 - sin^ (26) 



and exhibits a typical oscillatory behaviour. The formula ( p6D was first introduced by 
Pontecorvo |jl2[ in analogy to what happens in the K'^ — K'^ system and, as we have 



shown, can be derived for neutrinos with both Dirac and Majorana mass terms. From 
(p6|) we emphasize that both m^) and tum must be non- vanishing for neutrino- antineutrino 
oscillations to occur. 



3 Neutrino-antineutrino oscillations in matter 



When neutrinos propagate in a medium, their self-energy, acquired through the interaction 
with the particles in the medium, must be taken into account in writing the equations of 
motion; substantially, this is done with the replacements of the mass terms with the self- 
energy terms, as explained for example in For non- magnetized media (the extension 
to media with a magnetic field is straightforward following [0), at first order in the Fermi 
coupling constant Gp, the self-energies are given by [|13[, 



K 





-bL 




75 



2 

1 + 75 



(27) 
(28) 

(29) 
(30) 
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where is the medium 4-velocity (we will consider the medium rest frame, = (1, 0)) 
and 

-bL = V2Gf (n, - l-N^ (31) 



for the electron flavour, while 

- 6l = ^ iV„ (32) 

for the fi and r flavours, N^. and A'"„ being the electron and neutron number density of 
the medium respectively. The equations of motions are: 

{/c - m+) + {/: - m+ + bi /i) i^L + (/^ - rn+) + 

+ {/i - m+ - bL A)^R = (33) 

- (/; - m_) ur - (/; - m_ + fei jh) ul + {J^ - rnS) + 

+ (/; - m_ - 6i ;^)z/g = (34) 

and, explicitely, in the chiral Weyl base. 



(jj — cr • k — m+ I \ Vl 



+ 



+ f -"^^ + = (35) 



— a ■ k — m+ j I 

m_ - (a; + + (T ■ k) W ///J 

(a; — cr ■ k) m_ / I 

a; — a ■ k — m_ ) \ vS. 



+ 

+ 1 , """'^"■M 1 = (36) 



Introducing now, as above, the helicity eigenstates, in the (well verified) approximation 
that the interaction with the medium does not change neutrino helicity, the equations 
(§1), (§g) take the form 

(a; + A A;) ni — m+ — b^nx = (37) 

— m+ni + {uj — Xk)n2 = (38) 

{u + Xk)ni — m_ n2 — biUi = (39) 

-m_ni + (to - Xk)n2 = (40) 

With simple algebra we can write these equations as 

- - ml) 712 = bL — (uj - Xk) n2 (41) 

(u"^ - k"^ - m^_) n2 = bL — (uj - Xk) n2 (42) 
^ ^ m+ 

u — Xk 

rti = n2 (43) 

m+ 

oj — Xk ^ , ^ 

ni = n2 (44) 
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These equations admit non trivial solutions only if 

(,^2 _ ^2 _ ^^2 _ p _ ^2^^ = bl{u - Xkf (45) 



Eq. ( p5D is the dispersion relation for Dirac-Majorana neutrinos propagating in matter. 

Let us now focus on A = — 1 case; in the ultrarelativistic limit, the dispersion relation 
(^) can be viewed as the eigenvalue equation for the "effective hamiltonian" 

in the base ( | . 

V "1 / 

The time evolution of this states is then obtained by the Schodinger equation 

= (47) 



dt \ ni \ n 



For the sake of simplicity, we will consider a constant density medium; then the hamilto- 
nian in p7|) can be diagonalized by a matrix 

^ ^ / cos e sin e \ 

\ - Sm Ora COS Orn ) 

(neglecting CP violating effects) depending on a constant mixing angle between rii and 
ni. The eigenstates are then 



nim \ ^ ( cosOm - sin 9m \( rii 
I V sin 0m cos Om j \ ni 



(49) 

and the matrix U'^ -f^e// U becomes diagonal when the mixing angle is given by 



tan2^„ = -^^I^ (50) 
niD ruM 



The time evolution of the physical matter eigenstates is 



nvm{t)> = e-*"-* |r^l„^(0) > (51) 
rh^nit) > = e-'^'\rHmiO) > (52) 



and the eigenvalues are: 
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From the time dependence of the weak-interaction eigenstate pl-, 

6-^"^+*)) \ y{Q) > + 

(54) 

we can finally obtain the expression for the survival probability for an initial produced 
P{vl ^l) = 1 - cos^ 29m sin^ ^ t = 

We immediately note, from (|55|), a completely different behaviour of neutrino-antineutrino 
oscillations in matter with respect to the case of flavour oscillations. In fact, matter 
neutrino-antineutrino oscillations have both an amplitude and a period (and then an os- 
cillation lenght) sma//er than those corresponding to vacuum oscillations. In other words, 
opposite to what happens for flavour oscillations, neutrino-antineutrino oscillations in 
matter are always suppressed with respect to the vacuum case; furthermore, no resonance 
condition (analogous to the MSW effect 0) can occur. 



(( 



^-iuj-t _|_ ^-iuj+t 



COS 29r, 



-i ui^ t 



+ sin 29m (e 
w{0) > 



•i U)- t 



4 Limits on the Majorana neutrino mass 

The fact that neutrino-antineutrino oscillations in matter are suppressed with respect to 
the vacuum makes this phenomenon pratically non-interesting for astrophysical scopes. 
In particular, there cannot be a matter solutions in terms of neutrino-antineutrino os- 
cillations to the solar neutrino problem [^]. This can be easily seen from the following 
arguments. The amplitude cos'^29m in (^) is lower for more energetic neutrinos (such as 
those produced by the Be^ and reactions in the sun |T^) than for less energetic ones 
(such as those from the p-p chain). So one predicts a smaller reduction for Be^ and B^ 
neutrinos than for p-p neutrinos. This is just the opposite of what is needed to explain 



the experimental results [1^, so the matter effects on neutrino-antineutrino oscillations 
cannot help to solve the solar neutrino problem. 

On the contrary, the analysis of the experimental data on solar neutrino [^] shows 
that a vacuum neutrino-=antineutrino oscillations solution is indeed possible with the 
following value of the product of Dirac and Majorana z/g masses: 

mo rriM ^ 6 x 10"" eV^ (56) 

Note that this solution is equivalent to the vacuum flavour oscillation solution with max- 
imal (flavour) mixing angle (sin^ 29 ~ 1) and Am^ = mf, — ~ 6 x 10"^^ eV^ 
ITT 



7 



However, the best probe for neutrino-antineutrino oscillations (in vacuum) is with 
laboratory experiments. 

We have analyzed the current data on disappearance laboratory experiments for giving 
limits onmDniM- The most stringent limit for the electron flavour comes from the exper- 
iment at the Krasnoyarsk reactor []18[; with the formula (|26|) for the survival probability 
we have obtained 



rrfj^ m^^j < 7.5 x 10"^ eV^ 



(90% C.L.) 



(57) 



For the muon flavour, experiments allow both low and large values for m^^ (see 
from the Particle Data Group analysis ||19| we deduce 



< 0.23 eV^ 



(58) 



or 

m'^j^m^i > 1500 (59) 

Given the form (|26|) for the survival probability, from the experiments one can only deduce 
limits on the product of the Dirac times the Majorana mass. However, if the neutrino 
Dirac mass terms are generated by the same Higgs doublet giving mass to the charged 
fermion (as happens in G.U.T.), it is natural to assume that rriD is of the same order of 
magnitude of the (Dirac) mass of the charged fermion associated to the given neutrino 
in the same supemultiplet (for example, in SO(IO) u^., f^, are linked to the up-quarks 
u,c,t). Assuming then m^^ ^ 0{MeV) and m^^ ^ 0(100 1000 Mel/), from (|8D, (|5|) 
one can deduce the following indicative limits on and Majorana masses: 

m^f^j < 10~^ eV (60) 

and 

m^/ < 10-^ eV (61) 

or 

m^./ > 10'*^ eV (62) 

Note that the limit for is smaller by several orders of magnitude than the ones 
obtained from neutrinoless double (3 decay (some eV) PD| . 



5 Conclusions 

In this paper we have analyzed the theory of Pontecorvo neutrino-antineutrino oscillations 
in vacuum and in matter, showing that they can occur only if neutrinos are Majorana 
particles with both Dirac and Majorana masses non- vanishing. With a fleld-theoretic 
approach we have rederived the Pontecorvo formula (eq. (0)) for vacuum oscillations 
and then generalized this formula to matter oscillations, flnding that neutrino-antineutrino 
transitions in matter are suppressed with respect to vacuum ones. 
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Furthermore, we have shown that for solar neutrinos, the matter effects are not in 
agreement with the experimental data, so that matter Pontecorvo oscillations does not 
explain the solar neutrino problem. However, a vacuum solution to this problem in terms 
of neutrino-antineutrino transitions indeed exists, and we give the corresponding range of 
the product of Pe Dirac and Major ana mass. 

The existing limits from laboratory experiments for this product for electron and muon 
flavours have also been analyzed, leading to limits on Majorana mass (in the framework 
of Grand Unified Theories) several orders of magnitude smaller than the one found from 
the experiments on neutrinoless double beta decay. 
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